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Specific heat data for the quasi one dimensional quantum 
magnet Copper Benzoate (Cu(C6DsCOO)2 • 3D2O) is ana- 
lyzed in the framework of an effective low-energy description 
in terms of a Sine-Gordon theory. 

PACS: 74.65. +n, 75.10. Jm, 75.25,+z 



I. INTRODUCTION 

Quasi one dimensional quantum magnets have been a 
focus of intense theoretical and experimental interest for 
a long time. Most of the work is based on and centers 
around the spin- 1/2 Heisenberg model 01. 

Starting with Bethe's seminal work B] a host of ex- 
act results have been obtained for ground state prop- 
erties ||, magnetic susceptibility Q, thermodynamics 
i[, excitation spectrum |7],|| and correlation functions 
12| . From the point of view of standard spin wave the- 
ory which is highly successful for "three dimensional" 
materials, the findings of the these investigations were 
rather unusual. Over the last thirty years a number 
of anisotropic materials have been found that constitute 
excellent realizations of the one dimensional Heisenberg 
model ||, e.g. KCuF 3 , Sr 2 Cu0 3 , Cs 2 CuCl 4 or CuPzN, 
and many theoretical predictions have been confirmed ex- 
perimentally. One main focus of attention was the spec- 
trum, which comprises of an incoherent (two particle) 
scattering continuum of elementary excitations, the so- 
called spinons 0. These can be visualized in terms of 
ferromagnetic "domain walls" and are strikingly differ- 
ent from the usual spin waves. In particular spinons are 
believed to have fractional (semionic) exclusion statis- 
tics The low-energy effective theory of the spin-1/2 
chain is simply a free massless boson p|Jl5|Jlq| 



C 



(1) 



It has been known for some time that Copper Benzoate 
is another realization of a quasi-lD S = 1/2 Heisenberg 
antiferromagnet |i7| . However, its response to a mag- 
netic field has been found to be unusual Jig] : structural 
anisotropy leads to to generation of small staggered fields 
in and perpendicular to the direction of the applied field. 
Early specific heat measurements [ fl9| showed behaviour 
incompatible with theoretical results for a simple Heisen- 
berg chain. 

In a series of recent experiments |2C].2l|] the behaviour 
of Copper Benzoate in a magnetic field was investigated 
in great detail. Neutron scattering experiments |HJ es- 
tablished the existence of field-dependent incommensu- 



rate low energy modes. The incommensurability was 
found to be consistent with the one predicted by the 
exact solution of the Heisenberg model in a magnetic 
field. However, the system exhibited an unexpected ex- 
citation gap induced by the applied field. As no evidence 
for ordering was found in the experiments down to tem- 
peratures of 0.1 — 0.3isT the interchain coupling in Cop- 
per Benzoate supposedly is very small (the exchange is 
approximately 18K). We therefore will neglect it in the 
present work. 

In p^] it was proposed that Copper Benzoate is de- 
scribed by the Hamiltonian 

ffcuB = ■ $+1 - E»bHS? + nuhi-iySf , (2) 

i 

where H ^S> h, s is the effective Lande g- factor and 
,/ = 1.57 meV (g0[. Here the induced staggered field 
h is a function of the known (staggered) g-tensor |T§| ] 
and the Dzyaloshinskii-Moriya (DM) interaction in Cop- 
per Benzoate, for which unfortunately only scant infor- 
mation is available. If direction and magnitude of the 
DM interaction are given, h is calculated as follows [ p2| . 
The <?-tensor in the a" , b, c" basis (these denote the three 
principal axes of the exchange interaction p3[]) is given 

by @ 



9 = 



2.115 
±0.0190 
0.0906 



±0.0190 0.0906 
2.059 ±0.0495 
±0.0495 2.316 



(3) 



The ± correspond to the two inequivalent sites and in- 
dicate that application of a uniform field induces a stag- 
gered one. The corresponding contribution to the Hami- 
tonian is 



Hmagn — 9ab M-B H a S . 

The staggered DM interaction is 

Kdm = ]T(-1) J D ■ (Sj x S j+1 ) 



(4) 



(5) 



where \D\ <§; J and the direction of D is thought to be 
close to the a" axis. The DM interaction is eliminated by 
a rotation in spin space around the Z3-axis by an angle 
a = ±(arctanZ?/J)/2 on even/odd sites. This induces a 
very small exchange anisotropy which is negligible, and 
a staggered field 

Hoa = ~Hb(H x D/J) ■ ]T(-lM , (6) 
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where we used that ID|/J <C 1. Combining the two con- 
tributions (H) and dj) we obtain the total induced stag- 
gered field. For example, a uniform field applied along 
the a" axis induces a staggered field h along the b axis 
of magnitude (0.019 - 2.115D c /2 J)H. 

The low energy effective theory of @j is obtained 
by abelian bosonization and is given by a Sine-Gordon 
model with Lagrangian density [E2[ 



2M sinmr£/2 ; n = 1, . . . , [l/£] 



(11) 



C=l(d^) 2 + \(h)cos(l3 6) 



(7) 



Here is the dual field and the coupling (3 depends on the 
value of the applied uniform field. The coefficient A(/i) 
can at present not be calculated exactly. The reason 
is that the amplitudes of the bosonized expressions of 
lattice spin operators for H > are not known (in the 
absence of a magnetic field they have been determined 
very recently in |l^,^,^5|). For later convenience we 
define the coupling 



ln-f3 2 



(8) 



The Sine-Gordon theory (^|), for all its apparent simplic- 
ity, has fascinated physicists for decades. It is of interest 
as an integrable classical nonlinear differential equation 
featuring soliton solutions. On the quantum level it has 
been one of the cornerstones of nonperturbative quantum 
field theory with many exciting features such as quantum 
solitons, topological charge or regularization dependence 
in the nonperturbative regime [ p6[ . Most importantly 
the quantum Sine-Gordon model is exactly solvable |27| | 
and many physically important quantities can be calcu- 
lated. In particular, the spectrum is known to consist of 
a soliton-antisoliton doublet of mass M and their bound 
states which are called "breathers" . 

The soliton mass gap M can in principle be calculated 
exactly pi| in terms of (3 and A for a given short-distance 
normalization of correlation functions. However, A is 
known only for the case of vanishing uniform field H = 
(see above). A simple analysis based on the results of 
p8| , ^6|j25| yields for this case (i.e. one takes to staggered 
field into account, but bosonizes at H = 0) 



M/J 



2r(i/6) 

v^r(2/3) 

^2/3 

J 



£(3/4) 
T(l/4) V 4(2^) 3 / 2 



2/3 



2/3 



1 



(9) 



where we have neglected logarithmic corrections. This 
is in good agreement with the numerical analysis of the 
lattice Hamiltonian (0) for H = 0, which gives |25| 



M/J sa 1.85 
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(10) 



In addition to soliton and antisoliton there are [|] (here 
[ ] is the integer part) breathers with masses 



The mass spectrum as a function of magnetic field for 
Copper Benzoate was explicitly determined in ]29|] . 

Exact predictions of the low-energy effective theory (Q) 
for the spectrum p2| and the dynamical structure factor 
p9| were found to be consistent with Neutron scatter- 
ing experiments for applied fields along the b-axis. It 
is interesting to note that the Sine-Gordon solitons and 
breathers are fundamentally different from the spinons of 
the spin- 1/2 chain. 

In pi] ] precise measurements of the low-temperature 
specific heat were presented and analyzed in terms of sev- 
eral noninteracting one-dimensional bosons of the same 
mass. On the other hand, the spectrum of the Sine- 
Gordon model in the relevant region of couplings features 
five interacting modes with different masses fl29| (soliton, 
antisoliton and three breathers). 

In the present work we analyze the specific heat data of 
|ilf in the framework of the Sine-Gordon theory. A very 
important input in the low-energy effective Lagrangian 
(0) are the values of the coupling f3 and the spin velocity 
v s . In a "minimal" model (MM) they are calculated from 
the exact Bethe Ansatz solution of the Heisenberg XXX 
chain in an applied magnetic field H |l5| ] . In Appendix A 
we summarize the corresponding relevant Bethe Ansatz 
results. 
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FIG. 1. Coupling constant 2 /2tv in the MM as a function 
of the applied magnetic field H || b . 
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FIG. 2. Spin velocity in the MM as ; 
field H || 6 . 



function of the applied 



This procedure appears to be reasonable as long as 
the induced field h is very small so that its effects on 
and v s are negligible. The results are shown in Figs |l| 
and |2| respectively. At very high fields g^BH ~ 23 we 
approach the incommensurate-commensurate transition 
to the saturated ferromagnetic state |30| and the spin 
velocity thus tends to zero. 

An alternative scenario is to consider the spin velocity 
and/or the coupling as a phenomenological parame- 
ters in the Sine-Gordon theory §7§. The rationale behind 
such an approach is that the known presence of the DM 
interaction as well as possible exchange anisotropics will 
lead to deviations in the values of these quantities as 
compared to the MM predictions. 

A simple calculation shows that the effect of the DM 
interaction on f3 and v s is negligible. Adding an exchange 
anisotropy 



-ffanis = — JA ^ SjS* 



'j+i 



,0< A< 1 



(12) 



to the Hamiltonian (|^) firstly induces a change in the 
spin velocity entering the effective Lagrangian (0) and 
secondly generates the second harmonic of the SG inter- 
action, i.e. the effective low-energy theory becomes J3l| 



Acos(/3 8) +^cos(2/3 9) 



(13) 



Here we have assumed that A is much smaller than the 
magnetic energy scale g^-sH ■ The coupling /i mainly de- 
pends on A (the second harmonic also gets generated at 
1-loop level by the cos /38 interaction) . In the regime of 
couplings f3 we are interested in, the second harmonic is 
a relevant operator (in the RG sense), although it is of 
course much less relevant than cos/? ( I>. This means that 
we can safely neglect the second harmonic, unless A <C /i 
! We will return to this point below. Physically the effect 
of (na) is the following: if no uniform field is applied, the 



system remains critical. The spin velocity and the crit- 
ical exponents are changed slightly. If a uniform field is 
applied perpendicular to the direction of the anisotropy, 
a spectral gap forms, even if no staggered field is gener- 
ated. 



II. SINE-GORDON THERMODYNAMICS 

The thermodynamics of the Sine-Gordon model is 
most efficiently studied 32 via the recently developed 
Thermal Bethe Ansatz approach [1331] , which circumvents 



problems associated with solving the infinite number of 
coupled nonlinear integral equations that emerge in the 
standard approach based on the string hypothesis |34|] 
(note that the coupling constant (3 in our problem is a 
continously varying quantity and no truncation to a finite 
number of coupled equations is possible). It was shown 
in |32j that the free energy of the Sine-Gordon model can 
be expressed in terms of the solution of a single nonlinear 
integral equation for the complex quantity e{9) (we set 
the spin velocity to 1 for simplicity) 

e{9) = -iMpsinh(9 + ir]') 

dd'G (9 - 6') In (1 + exp \-e{9')}) 

+ / d6'G (9-6' + 2irf) In (1 + exp [-e(0')D , (14) 



(15) 



where f3 = (feT) 1 , M is the soliton mass and 
dui cos(2oj9/tt) sinh(cj(£ — 1)) 



G (9) 



sinh(w^) cosh(aj) 



The free energy density is given by 
f((3) = -—Qm / d9 sinh(6> + irf) In 

P7T J-oo 



1 + e 



-6(6) 



(16) 



As we are interested in the attractive regime 7 < 4-7T 
we have 



< 77' < tt£/2 



(17) 



Note that the free energy does not depend on the value 
of r/' as long as it is chosen in the interval (|l7]). The set 
( |l4| ) of two coupled nonlinear integral equations is solved 
by iteration. For j3 — * 00 the first iterations can be cal- 
culated analytically and the corresponding contributions 
to the free energy are seen to be of the form 



fW 



2M 
Mi 



LXJ / 

E- 

71=1 



-Ki(nm(3) 



KM/3) 



(18) 



where Mi — 271/8^^ is the mass of the first breather 
and K\ is a modified Bessel function. The first term is the 
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contribution of soliton-antisoliton scattering states to the 
free energy, whereas the second term is the contribution 
of the first breather. Both terms have the form character- 
istic of massive relativistic bosons. The contributions of 
the heavier breathers are found in higher orders of the it- 
erative procedure employed in solving (|lj). The specific 
heat is obtained from the free energy 



C — T 



g T 2 



(19) 



At low temperatures it is found to be of the form 



C 



[i/«] 

E 

a=l 



27TU, 



k B T 



k B T 



3/2 



eM~M a /k B T) , 



(20) 



where M a are given by In order to compare the- 

oretical predictions based on the SGM with the specific 
heat data of [Ell we need the free energy at "intermedi- 
ate" temperatures and thus have to resort to a numerical 
solution of ( |l4|) by iteration. 



III. SPECIFIC HEAT IN COPPER BENZOATE 

Let us now investigate the question how well the theo- 
retical predictions based on an effective Sine-Gordon the- 
ory agree with the specific heat data of pjf . As was 
pointed out in pl| , at very low temperatures a nuclear 
contribution to the specific heat is present. In the follow- 
ing analysis we neglect this contribution but note, that 
by taking it into account we can achieve excellent agree- 
ment of the SG results with the data at low temperatues. 

We now analyze the specific heat data of [^l] as follows: 
we calculate the specific heat in the framework of the MM 
using the soliton gap as a free parameter, which is then 
fixed by fitting the calculated specific heat to the data. 
This procedure yields the dependence of the gap on the 
applied field M(H), which has to be consistent with (|1C) 
and the dependence of h and H, which follows from (M) 
and (^). In order to keep things simple we ignore the 
logarithmic correction and independence of (3 in ( |l0| ) so 
that 



M(H) w cH 2/3 . 



(21) 



Here the coefficient c depends on the orientation of the 
applied field and the direction and magnitude of the DM 
interaction as explained above. In order to calculate c, 
we would need to know the precise magnitude as well 
as orientation of the DM interaction as is clear from @ . 
Unfortunately this information is presently not available. 
From considerations based on the crystal structure D is 
expected to lie in the a" — c" plane and and is thought 
to be roughly of magnitude D/ J ~ .1. 



Reversing the logic [E2| , if we determine the coefficient 
c in (|2l|) for all three independent orientations of the 
applied field by fitting the SGM predictions to the data, 
we can calculate what the direction and magnitude of D 
has to be in order to reproduce these results. Below we 
follow this line of argument. 



A. Magnetic Field along c" axis 

For magnetic fields along the c" axis we find excellent 
agreement of the data with the "minimal" model dis- 
cussed above. This is shown for some values of H in Figs 
| and |. 




FIG. 3. Specific heat as a function of temperature for fields 
of H=1T and H=2.5T applied along the c" axis . 

As explained above, the presence of an exchange 
anisotropy would change the spin velocity. Assuming v s 
to be eight percent smaller than in the MM we still ob- 
tain good agreement with the data (note that the soliton 
mass is of course changed as well) as is shown in Fig. ^ 
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FIG. 4. Specific heat as a function of temperature for fields FIG. 6. Specific heat as a function of temperature for fields 

of H=3.5T and H=7T applied along the c" axis . applied along the b axis . 




FIG. 5. Specific heat as a function of temperature for fields 
of H=3.5T and H=7T applied along the c" axis. The spin 
velocity is taken to be eight percent smaller than in the MM 



In order to check which scenario for v s is in better 
agreement with experiment the values for the soliton 
masses obtained by fitting to the data have to be com- 
pared with (|l0|). 



B. Magnetic Field along b axis 

For magnetic fields along the b-axis the agreement of 
the MM prediction with the data is less impressive. As 
is clear from Figs the MM systematically underesti- 
mates the measured specific heat in the temperature re- 
gion T w 0.8 — IK although there is still fair agreement 
of the MM with experiment. 
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FIG. 7. Specific heat as a function of temperature for fields 
applied along the b axis . 

A much improved fit to the data is obtained if the 
spin velocities in the effective SGM are changed by eight 
percent as compared to the MM. This is shown in Figs ||- 




FIG. 8. Specific heat for fields applied along the b axis. 
The spin velocity is taken to be eight percent smaller than in 
the MM . 
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FIG. 9. Specific heat for fields applied along the b axis. 
The spin velocity is taken to be eight percent smaller than in 
the MM . 

In order to check the compatibility of the fitted values 
for the soliton gap Mat(H) with (^|) we plot M&t(H) as a 
function of H in Fig. [HI For simplicity we only consider 
the results calculated on the basis of the MM. We find 
good agreement for applied fields along the b and c" axes. 
The logarithmic corrections ( |To| ) to the gap may improve 
the agreement, but need (unavailable) information on the 
DM interaction as input. The ratio of mass gaps for fields 
applied along the b and c" axes is found to be 



A c »/A 6 re 1.43/0.65 = 2.2 



(22) 
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FIG. 10. Gap of the soliton for H || c" and H \\ b. The fits 
are to the simple scaling law M oc (H/J) 2 ^ 3 as described in 
the text . 



C. Magnetic Field along a" axis 

For fields applied along the a" axis it is impossible to 
obtain agreement of the MM predictions with the data. 
How can we understand this fact? The DM interaction is 
expected to lie in the a" — c" plane, so that we can write 



D 



D c e c n 



(23) 



For H — He a " the net induced staggered field is directed 
along the b axis and is of magnitude 



h = (0.019 - 2.1 15D c /2 J)i? 



(24) 



Clearly this would be very small if D c /J pa 0.02. We 
note that such a value of D c together with the gap ratio 
A c «/A 6 ps 2.2 implies that D a /J ps 0.12. Note that this 
is consistent with the expectation that D/Jm .1 and a 
direction of D close to the a" axis. 

In this case the coupling X(h) in would become 
very small and there would be a regime in which pertur- 
bations other than the staggered field would dominate 
the physics. For example, if an exchange anisotropy jl^ ) 
was present in Copper Benzoate, the low-energy effective 
theory would be given by ( |l3|) with /i ^S> A. As a first 
approximation we then can ignore the A cos (3® term and 
study the remaining (repulsive) SGM. In the framework 
of this scenario we obtain a rather reasonable fit to the 
data as is shown in Fig. The expected mass gap is dif- 
ficult to estimate, because as we already mentioned the 
coefficients in the bosonization formulas are known only 
in the absence of a magnectic field [^6| , ^5|j2^ 1 . A crude es- 
timate can be obtained by approximating the coefficents 
in the presence of a uniform field by the ones for H = 0. 
The gaps obtained by fitting the data are found to be 
consistent with a rather small anisotropy A/J^0.05 in 
this approximation. We note in passing that the zero- 
field specific heat found in pl[ 



C(T) = 0.68(l)i?fc B T/J 



(25) 



actually corresponds to an anisotropy of the type ( |l2| ) 
with A/ J ps 0.06. 

In order to work out a more quantitative theory for 
the a"-axis specific heat data the full two- frequency Sine- 
Gordon theory would need to be analyzed, which is possi- 
ble in a perturbative framework |35| . We hope to address 
this point in the future. The (non) existence of the mech- 
anism described above could in principle be checked by 
inelastic Neutron scattering with H \\ e a >>: if the physics 
is indeed dominated by interactions other than the in- 
duced staggered field, the spectrum will be very different 
from the one oberved in |pl| . In particular, if exchange 
anisotropy is the relevant mechanism and the effective 
low-energy theory is thus given by (U) with A ps 0, then 
no coherent one-particle excitations are present. The dy- 
namical structure factor at wave number tt is then dom- 
inated by an incoherent soliton-antisoliton continuum. 
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FIG. 11. Specific heat for H \\ a". The theoretical curve is 
obtained by assuming an exchange anisotropy as explained in 
the text . 



IV. CONCLUSIONS 

We have analyzed specific heat data for Copper Bcn- 
zoate in the framework of a Sine-Gordon low-energy ef- 
fective theory. For uniform magnetic fields applied along 
the b and c" axes we find good agreement of the theory 
with the specific heat data. The a" axis data cannot be 
understood by the same theory that applies for the b and 
c" axes. We argue that the staggered field induced by 
the DM interaction essentially cancels the field induced 
by the staggered g tensor for H || e a « so that a new mech- 
anism is responsible for generating the gap. We propose 
that exchange anisotropy might be responsible. 
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APPENDIX A: 



THE HEISENBERG CHAIN 
IN A FIELD 



We summarize some relevant results (for derivations see 
|^5| ) for the anisotropic Heisenberg model in a magnetic 
field 



Hxxz = 



S n+1 + (A - l)S*S* n+1 - 2hS z 



(26) 



where — 1 < A = cos 7 < 1 and J > 0. The customary 
form of the Hamiltonian is obtained by performing the 
unitary transformation 



5 ,, 



6V, 



(27) 



At low energies (|26| ) is described by a free massless boson 
(|l|) compactified on a ring of radius R i.e. $ and $ + 2irR 
are identified. The dual field 9 fulfils 6 = 6 + 1/R (see 
e.g. j36|). The following bosonization rules can be derived 
along the lines of e.g. |37): S n — > [J(x)+n(x)], x = nao, 
where ao is the lattice spacing and 

r- = ^d x $(x) , 
p 

(-1)"J+ =iA(6x)exp(if3Q(x)) a d x G(x) 

exp (— i/3Q(x)) sin (~^®( x ) ~ ^x^j , 



n x {x) = A ± cos(f30{x)), n y {x) = A ± sin(/36(cc)), 



n z (x) = -(-l) n A z sin(^$(x) - 26x). 



(28) 



Here (3 = 2nR and the coefficients Aj^ z are known only 
in the absence of a magnetic field The standard 

structure in terms of uniform and staggered magneti- 
zation operators is obtained by performing the unitary 
transformation (f27j). We note that the often neglected 
first term in J + is actually more important than the sec- 
ond: as a matter of fact, for H = it yields the leading 
contribution to transverse correlations at wave number 

7T of (||) fUH. F ° r & = ^ Le ' the invariant 
chain in zero field, the second term is simply the sum of 
left and right SU(2) currents. The first contribution to 
J + corresponds to a particle-hole excitation with spin 1 
relative to the ground state (see (II. 5. 3) and (XVIII. 1.16) 
of jl5]]) and for H = can be derived by carefully tak- 
ing the continuum limit of the Jordan- Wigner lattice 
fermions p9| , ^6| . For H 7^ we expect A(Sx) on general 
grounds to be of the form const x cos(2Sx). Equations 
(|28| ) are used to derive the continuum form of the pertur- 
bation (|l^) (note that both smooth and staggered com- 
ponents of the spin operators contribute for (3 — V2tt). 

The constant /? and spin velocity v s are determined by 
A, J and h of the lattice model as follows. The dressed 
energy e, momentum p and "charge" Z of an elementary 
spinon are given in terms of the solutions of the linear 
integral equations 



r(A) 



^K{\~n) e{jx)=2h- 



4 J sin 7 
cosh2A -I- cos 7' 



A 



p{\) = 2tt dn p(n) , 
Jo 

; _ A27r *(A-M) ^- 27r [ cosh 2A + cos7] ' 

rA 



P(A) 
Z(X) 



2^ 

dfj, 
2^ 



2 sin 7 



K(X - /x) Z(jm) = 1 , 



(29) 
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where K(X) = 2 sin27/(cosh A — COS27). Here A is the 
rapidity corresponding to the Fermi momentum and is 
fixed by the condition 



[5] 



e{±A) = 



(30) 



The spin velocity is then given by the derivative of the 
spinon energy with respect to the momentum at the 
Fermi surface 



de(X) 



dp(X) 



de(X)/dX 



\=A 



2tt P {X) 



\=A 



Finally, (3 and 5 are given by 



/3=-^,£=~-7r/ d[i p(fj,) 



(31) 



(32) 



In order to determine v s and (3 we solve ( p9| ) numeri- 
cally, which is easily done to very high precision as the 
equations are linear. The results are shown in Figs ^ and 
^|. Finally, we note that correlation functions at small fi- 
nite temperatures can be calculated like in (h]] (see also 
p0|). We only must remember to shift the momentum 
by ±2(5 away from n for the longitudinal correlation func- 
tion and use the correlation exponent as calculated above 
from the Bethe Ansatz. For example the transverse dy- 
namical susceptibility at small momentum (which cor- 
responds to momentum n in the customary form of the 
Heisenberg Hamiltonian, which is related to ( p6f ) by the 
unitary transformation (p7|)) is given by 



,uj + v s q 



v s q] 



Sn 



4ttT 



4ttT 

2 
1 - — ) 



2 
4ir 



(33) 



where B(x, y) is the Beta function and q is close to zero. 
The longitudinal susceptibility is dominated by the gap- 
less modes at 7r ± 28. It is the sum of two terms 



<j=± 



' 2 ' 



2/3 2 1 47rT ' J ai) 



.U! + V s Qo 



(34) 



where Q a — q — it + a2S. 
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